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Abstract 

The wave function of a light pscudoscalar meson is considered and nonpertur- 
bative corrections as signaled by perturbation theory are calculated. Two schemes 
are used, the massive gluon and the running coupling scheme. Both indicate the pres- 
ence of leading power corrections of 0(6^), whose exponentiation leads to a Gaussian 
dependence of the wave function on the impact parameter b. The dependence of 
this correction on the light cone energy fractions of the quark and the antiquark is 
discussed and compared with other models for the meson. 
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1 Introduction 



The challenging task of explaining experimental data for exclusive processes, such as pion and 
nucleon electromagnetic form factors and elastic scattering, is tightly connected with our under- 
standing of the bound states of light hadrons. In the asymptotic limit of parametrically large 
momentum transfer (Q) there is a clear theoretical picture that comes from merging the parton 
model with perturbative QCD (pQCD). The basic tenet of this picture is that a hadronic bound 
state is a superposition of virtual states with a definite number of constituent partons. In the frame 
where all participating hadrons are fast moving, factorization guarantees that the dependence 
of an exclusive observable enters only through the underlying process of the elastic scattering of 
the partonic constituents. This has lead to the, by now classic, dimensional (or quark) counting 
rules Q. However, the fact that the scaling with as predicted by the counting rules is actually 
observed by experiments at high momentum transfers generated new questions. Counting rules 
assume that the dominant configurations are those in which none of the constituents carries van- 
ishingly small longitudinal momentum. But such configurations are present in the hadronic wave 
function and would lead to scaling with that is less steeply falling as ^ oo. Underestimating 
the contributions from these end-point regions was the main criticism against the parton model 
picture [^. 

The next major step forward came with the work of Sterman and collaborators |^ Q who im- 
plemented the summation of leading and next to leading logarithmic radiative corrections through 
the introduction of Sudakov factors into the factorized expressions for exclusive processes. The 
generic form of such factors is exp [—clnQ'^ ln(ln Q^/ ln(l/6^))] where b is the transverse distance 
between constituents. It has the effect of suppressing configurations of constituents that are sepa- 
rated by large distance 6, which would be the case in the end-point regions mentioned previously. 
We note that the Sudakov factors can be considered as the perturbative tail of the hadronic wave 
function, i.e. the region where the bound state properties can be reconstructed purely from pQCD. 
Sudakov improved perturbation theory, implemented with a model for the hadron wave function 
at low momentum scale and a reasonable prescription on how to freeze the coupling at large scales, 
does give numerically stable results with the smallest possible number of phenomenological pa- 
rameters [^, H, ^. This approach sets the benchmark but it is not the final answer. The main 
problem now becomes to estimate the sensitivity to soft contributions of the Sudakov improved 
perturbation theory for moderate Q^, that is in the range 3GeV'^ < < 40Gev^. To this end we 
need a framework for studying nonperturbative corrections to hadronic wave functions. This is the 
objective of this paper. 

Since our approach to the problem is from the high end, we consider the two quark 
component of a light pseudoscalar meson wave function. This is the dominant configuration in the 
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asymptotic limit. The specific exclusive process in which the pion participates will be of no concern 
here. Our goal is to calculate the structure of the nonperturbative corrections as they are signaled 
by perturbation theory itself. These will occur in the form of power corrections in b in the Fourier 
transformed wave function. 

Recently, some progress has been made in understanding the power corrections by consider- 
ing classes of Feynman diagrams that give rise to a factorial divergence of the perturbation series in 
large orders. For the case of QCD the power corrections arise from the small momentum region of 
the loop integrations and are associated with the so called infrared renormalons. For a review and 
related references, see [0, Thus, one may get an indication of the type of power corrections by 
looking at the infrared sensitivity of Feynman diagrams. Most investigations involving renormalons 
are in the context of two calculational schemes. The first is the massive gluon and the second is 
the running coupling scheme. In the first method, which we apply to the case of the radiative 
corrections to the meson wave function, a gluon mass is introduced via a dispersive parameter A 
1^] that acts as an infrared regulator. This procedure is certainly consistent at the one loop level. 
Since the power corrections, in the renormalon approach arise from the infrared sensitive regions, 
our interest will be in the structure of the nonanalytic terms in which, to the order we work in, 
turn out to be InA^ and ft^A^lnA^. It is known Q that such nonanalytic terms come only from 
the pinch singular points in the loop momentum integration with InA^ being associated with the 
usual logarithmic enhancements of the perturbative series and 6^A^lnA^ with the leading power 
corrections. In the second scheme, power corrections are computed using a one loop improved 
perturbation theory where infrared effects are introduced through the running coupling as{k'^). 
Then the Sudakov exponent is Borel transformed and its singularity structure in the Borel plane is 
studied. The power corrections in this case are found to be proportional to This technique 

is the same one used for the study of power corrections to inclusive processes and has been widely 
discussed in recent years, see for example Ref. [p^ . 

In section 2 we calculate the leading power corrections and discuss the correspondence 
between the expressions derived in the two schemes. Exponentiation of leading power corrections 
naturally lead to the Gaussian dependence of the meson wave function on the impact parameter 
b. In section 3 we discuss our results and compare them with previous publications where such 
Gaussian factors for the meson wave function have been advocated. Here, of interest will be the 
X and InQ^ dependence of the Gaussian factor itself, where x is the light cone fraction of the 
momentum of the quark in the meson. Finally we summarize our conclusions in the last section. 

2 Radiative corrections to the meson wave function 
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2.1 Definitions 

The Bethe-Salpeter two quark wave function of a fast moving pseudoscalar meson M is defined by 
the following matrix element at renormalization (and factorization) scale ^ 

x{k,p,^^) = ^J ^e*M0|r(g-(0)^75g(y))|M(p)). (2.1) 

In the frame where the meson is moving fast we define the light cone vector along the direction 
of motion = (l/-v/2)(l, 0, 0, 1) and the parity reflected direction vector yf^ = (l/\/2)(l, 0, 0, — 1), 
fj, = 0,1, 2, 3, normalized as v ■ v = 1. Then the light cone momenta of the meson and the quark 
constituent are defined as 

= p+jM ^ ^y.]^^ ^p+^ ^- ^ ^ . A;^ = k. (2.2) 

The large scale is p'^ = Ql\f2, where Q is the momentum transfer of the hard exclusive process 
in which the meson participates. Although we do not need to specify the process, we must always 
keep in mind that our discussion will always be in the context of a hard process. The object of 
interest in this work is the wave function in impact parameter space given by 



V{x,b,p,n) = J (fke'^-^ J dk-X{k,p,n) 
1 fdy 



Nc J 27r 



y (0|r(g-(0) /;75g(0,y^,b))|M(p)). (2.3) 



The quark distribution amplitude, which enters the leading order perturbative expressions, is 

(/>(x,/i)= ''d^k / dk- X{k,p,ii) =V^^\x,b = 0,p,fi). (2.4) 



Radiative corrections to V^^^ coming from the infrared region exponentiate. Their summation leads 



to the Sudakov suppression factor |11] 



Vix,b,p,fi) =expi-S^fix,b,Q,fi)) V^''\x,b = 0,p,ii), (2.5) 

with o o V 

2Cf / InQVA^CD \ 11 2 , , 

where NL denotes next to leading logarithmic corrections, which are also known. Note that the 
above perturbative answer is really defined in the region 1/6^ ^ Aqqj-,, i. e. the evolution of the 
wave function can be reliably constructed only for small size bound states. 

The Sudakov exponent shown in Eq. ( p. 61) is obtained from the one loop radiative corrections 
to 'P*^'^\ There exist various methods of calculating these corrections. In this paper we will mostly 
follow the approach of Botts and Sterman ^] . In this approach the fermions are taken massless and 
the axial gauge is used. The advantage of using the axial gauge in the derivation of an evolution 
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equation for the wave function is that the diagrams involved in the calculation of the kernel have 
simpler structure than in a covariant gauge. Specifically, diagrams in which gluons couple the kernel 
with the 'inside' of the wave function thus generating non-2PI connections in the constituent quark 
channel (non-ladder gluons) are suppressed, Fig. 1. In the covariant gauge such diagrams must 
be included and further reduced through the use of Ward identities, just like in the treatment of 
parton distribution functions in D.I.S. 



Fig.l Class of diagrams that contain 'non-ladder' gluons (/). These are suppressed in the axial 
gauge. 

In the axial gauge there are gauge (n) dependent contributions to S^f that enter in the next 
to leading logarithmic corrections in Eq. ( |2.6| ). Such gauge dependent pieces cancel against gauge 
dependent terms in the radiative corrections to the hard scattering. Since we do not analyze the 
hard subprocess here, our results will be at the level of leading logarithmic and power corrections. 

2.2 Meson wave function in the massive gluon scheme 

As was mentioned in the introduction, we are interested in the nonanalytic dependence of the 
radiative corrections to the meson wave function on the infrared cutoff. The Landau equations 
guarantee that infrared sensitivity arises from those momentum regions where the internal lines in 
a Feynman graph approach their mass shell. To regularize mass-shell divergences we introduce for 
the gluon a (dispersive) mass parameter A and calculate the one loop radiative corrections to 
in this massive gluon scheme. This is the main difference with the treament found in Ref. |^ 
and also the fact that we are mainly interested in the terms that are nonanalytic in A and vanish 
as powers of A in the exact mass-shell limit. The gluon polarization tensor N^^i^y is defined via the 
principal value prescription for the unphysical singularities n ■ q. Hence the gluon propagator reads: 



D^^uiq,n) = . N^^{q,n), 

q'^ — X'^ + le 



n ■ q m? 



N^,u{q,n) = -g^, + _ ^ n^^q,) - _ ^ q^q, , (2.7) 

with 7/ — > 0. At the end of this subsection we comment on the possibility of using the same parameter 
A to regularize the unphysical n ■ q singularities. The graphs to be calculated are depicted in Fig. 2 
and the Dirac structure of the meson-g-g vertex is X = —(1/4) /'75X, where X is the scalar wave 
function defined in Eq. (2.1). 
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(oi) (02) 
Fig. 2 One loop diagrams for the meson wave function V. 

Contribution to V from vertex correction, Fig. 2 graph(6): 
The vertex correction to the meson wave function is 



Vip) = i{4:TTas)CF J 



tr 



/'75 



1 



i)+ /k- /i + ie ]k- /i + ie 



-r 



N„ 



1 



(2.8) 

The impact parameter b acts as an ultraviolet regulator so the loop integral can be kept at = 4 
dimensions. Standard manipulations for the —g^u piece of the gluon propagator yield 



7^1 



as 



CfV 



(0) 
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+ 0(6^A^ln(6A)) +0(ln°(6A)). (2.c 



Without the gluon mass regulator, the ln(6^A^) term would turn into an IR divergence that should 
be dimensionally regularized. Here, as well as in the rest of the calculations the interest lies in the 
presence of the iP'}? ln(6^A^) nonanalytic term. 

The calculation involving the n-dependent pieces of the gluon propagator is somewhat 
subtler. Evaluation of the traces in Eq. (|2.8|) yields 

z(47ra,)Ci.p(°) {h + h) , 



with 



and 



d'^q jjj.q [n ■ {q — k) + V ■ nv ■ {q — k) — V ■ nv ■ q]n ■ ' 
(2^'' [{q - kf + ie] (g2 _ A2 + ie) [(n • qf + r?2] 

h 



(2.10) 
(2.11) 
(2.12) 



The two integrals Ii and I2 arise from the two pieces n^^j^ and riyq^ of the gluon propagator 
respectively. At this stage it is advantageous to fix the gauge by the choice 



(2.13) 



This gauge choice has the advantage of simplifying the corresponding integrals by removing all 



dependence on q and q, from the numerator of Eq. (2.11), i.e. the normal coordinates in the 
collinear limit for the gluon momentum. In this gauge Ii becomes 



h 



d^q 
(2^ 



[2n ■ q + {k ■ q) /(n ■ k) — 2n ■ k]n ■ q 
[{q - ky + ie] {q^ - A^ + ie) [{n ■ qY + if] ' 



(2.14) 
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The above integral is evaluated via the introduction of Schwinger parameters. The intermediate 
steps can be reconstructed by using the techniques and the integral relations presented in Ref. 1 12 1 
and the result is 

= TTyi t ^r'^' t [K^i^fB) + A^K^{^fB) \ , (2.15) 

with 

A = X2{-2 + X2 + xiX2)b'^{n ■ k)'^ , 



B = 



[1 - xi)x|(n • kf + (1 - X2)A2 + 



(2.16) 



Since we are interested in the limit A — > we will calculate leading and nonanalytic A-dependence 
of the Feynman integrals with the help of the Mellin transformation defined as 

/•OO 

M[F{t)]il3)= dtt-l^-^F{t). (2.17) 

JO 

It should be noted that the Mellin transform method has been used extensively in the study of the 
high energy behavior of Feynman diagrams. For a review of the method and for further references, 
see [|^]. The definition ( 2.17 ) yields the following correspondence between poles in /3 and (large) t 
dependence: 

r rff^o ln"~^ t 

«IW) = ^^^"F«) = ^^. (2.18) 

Nonanalytic terms in t (logarithmic) are generated by at least double poles in the Mellin image. 
We define the following dimensionless ratios 

t = i!i-|l!, , = 52(^.^)2^ (2.19) 
A^ 

and we note that the parameter P also acts as an IR regulator. Then, in Eqs. ( 2.16| ) the r] ^ 
limit for the gluon propagator can be taken. The Mellin transformation of the Bessel K functions 
can be found in Ref. |14| and the Mellin image of Eq. ( 2.15| ) turns out to be 

M[h] = j^iJo + Ji), (2.20) 

where 

MP) = 2^r(/3)s-^/2 /'dxixf/'(l-xi)-^/2 /'dx2xf(l-X2)-^K-/3(V5iTx2), (2.21) 

Jo Jo 

MP) = 2fT[P)s^l^'l"^ (\x,x-^'^^^/\l-x{)-^/^ 

Jo 

x/ dX2X^{l-X2)''^{-2 + 2x2-XiX2)Ki_p{./^X2). (2.22) 

Jo 

Inspecting the above integrals we see that for P > there are no singularities. Singularities arise 
only for P < and they are generated from the integration end points xi ^ and/or X2 0. The 
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extraction of such singularities is straightforward because the series expansions of the X-functions 
for small argument can be used. From Eq. ( p.l8 ) it follows that the relevant nonanalytic terms 
that we are interested in can arise from at least double poles at [3 = 0,-1 We examine these 
two points in turn. 

For the poles at /3 = we set /? = 5, 5 — > and obtain 

,1. 



Jo(/3 = 5) 



5>' 



Combining Eqs. ( |2.23 

-^1 1/3^0 



20|) and inverting the Mellin transformation we get 



(2.23) 



(47r) 



1 2 

2 {n-ky 



+ 2(1 -In 2) In 



(n • A;)^ 



+ R. 



(2.24) 



The residue R generically contains all terms that are analytic in A^. Similarly, the poles are (3 



are obtained by setting (3 



-1 + 5 and power expanding the i^- functions. The results are 
Jo(/3 = -1 + 5) = (s-independent) , (2.25) 



and 



-I + 5) = (s-independent) - ^ + ^(1 - ln2) + 0(1/(5) 



(2.26) 



Again, combining the above two equations with Eq. (2.20) and inverting the Mellin transformation 
we have 



(6— independent) 



1 



-h^\^ In^ 



A2 



(4^)2 
(6— independent) 



(47r) 



1 



In^ 



a2 ^ 



1 A^ 

-(l-ln2)62A2 In- -— 

2^ ^ (n-k)^ 



R 



+ -(1 - In 2) In 



A2 



X2Q2 



+ R. 



(2.27) 



In the last step we have used that {n-k)"^ = . By (6-independent) we denote all terms that are 
non-analytic in A^, 0{}? In A'^) in this case, and are independent of h. Analytic in A terms reside in 
R. We make this distinction because, as we will see later, the terms denoted by (6-independent) are 
of infrared origin and cancel against the self energy contribution. The expression for I2, Eq. (|2.12| ), 
is obtained from the above by simply — > (1 — x)^. The overall contribution to the vertex from 
the n(^gjy) piece of the gluon propagator is, Eq. ( 2.1C| ), 



as 

+0 



'ln2 + ln2 



A2 



x2Q2 



(1-X)2Q2 



+ (1 - In 2) In 



A2 



independent) + R . 



x(l -x)g2 

(2.28) 



^ It is straightforward to see that there are no other points in the range between -1 and that would give double 
poles, thus excluding any linear in A corrections arising from poles at /3 = —1/2. 
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Finally, the g^Qjy piece of the gluon propagator gives 6-independent leading contributions and hence 
it cancels against the corresponding piece of the self energy, as we will point out below. 

Contribution to V from self-energy correction. Fig. 2 graphs (oi) and (02): 
Due to UV divergences we employ dimensional regularization with D = 4 — 2e and renormalization 
scale /X. 



r, 



i{4TTas)CF ^J■ 



2e 



tr 



1 



1 



-7 



-7^ 



1 

"A2T 



(2.29) 



/c + ie ]k— /l + ie 

It is straightforward to find that, even before any gauge fixing, the self energy contribution is such 
that 

^(ai) + ^(a2) + V^h) (6^0)= Ruv , (2.30) 

where the residue Ruv contains terms of purely ultraviolet origin and it is due to the mismatch in 
the UV regularization of the two sets of diagrams. Note that the vertex diagram is UV regularized 
by a finite h whereas the self energy diagrams are dimensionally regularized. Had we taken the limit 
6^0 before calculating the vertex integrals, then this would have required to also dimensionally 
regularize the vertex and then in Eq. ( 2.30| ) Ruv would have been exactly zero. In any case, the 
net result of this analysis is that the self energy diagrams will cancel all 6-independent terms in the 
vertex, such as the ones shown in Eq. (2.28) as well as the vertex contributions from the q^Qu piece 
of the gluon propagator. 

The final result for the one loop radiative correction to the wave function in the massive 
gluon scheme is obtained by combining the partial results of Eqs. ( ^ , 2.28| ) and (2.30). Retaining 
only the nonanalytic in A terms up to 0{\^) we obtain: 



p(i) = ^Cf f Cih'X^ In + C2b'\^ ln2 X^)+R, 
47r V / 



(2.31) 



with 



Ci = -^\n{x{l-x)Q'), C2 = -. (2.32) 

The familiar Sudakov factor S^f ~ a^ln^Q^, calculated for fixed coupling a^, is analytic in A^ 
and it is contained in the residue R. In the coefficients Ci,C2 we have only kept the parts that 
are n-independent to leading logarithmic in Q order. Note that the n-independent leading contri- 
butions to Ci, C2 come only from Eq. (|2.28| ). The n-dependence enters through the combination 
ln(n • k)'^ + ln(n • {k — p))'^, and any variation of n will lead to change that is subleading in Q, i.e. 
it will be down by a ln(a;(l — x)Q^). Fixing the gauge as in Eq. ( |2.13| ) organizes conveniently the 
leading and subleading InQ^ contributions but the result that we quote in Eqs. ( 2.31 , 2.32 ) is the 
gauge independent leading piece. 

We now return to the question of regularizing infrared singularities in the axial gauge and 
within the massive gluon scheme. In the massive gluon propagator, Eq. ( p.Tj ), we used the parameter 
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A to regularize on-shell singularities whereas the polarization tensor N^y was defined through the 
principal value prescription and the regulator r/. If we treat the parameter A as a lagrangian mass, 
instead of a dispersive parameter (this is at the one loop level only) then the propagator will be 
modified by the addition of a fourth term of the form n^riy}? /[{n-k)'^ —X^in?]. A similar prescription 
has been considered in the leading A ^ limit in Ref. However, beyond the logarithmic in A 
level, this term yields 6-dependent 0{}? In'^ A^) contributions that cannot be interpreted as arising 
from some pinch singularity in the limit A ^ 0. Such artifacts of the regularization prescription are 
avoided once one keeps in mind that A is not a lagrangian gluon mass but just an on-shell regulator, 
as it appears in Eq. (^). 

2.3 The running as scheme and IR renormalons 

The running scheme for estimating nonperturbative corrections starts with the expression for 
the one loop radiative corrections to the meson wave function. The leading logarithmic term is 
well known and can be readily obtained from the expressions in Eqs. (^]^, 2 .291 ) after rationaliz- 



ing the fermion propagators and applying the collinear approximation to the numerator factors. 
Equivalenlty it can be computed by calculating emission from two, almost parallel, Wilson lines. All 
infrared divergences are dimensionally regularized and no A regulator need be introduced. Principal 
value prescription is used for the axial gauge gluon propagator which now reads 



q'^ + ie 



1 n2 

-g^lv H rll^q^,^ - q^qv 

n ■ q ^ ' (n ■ qy 



(2.33) 



After performing the collinear approximation to the numerator factors and then integrating over 
q~ by closing around its poles in the q~ complex plane we obtain the result Q 

m r<i^=Q^ d^ci / \ r^P'^ dn~^ 

V^'\x,b,p,^,) = ^CFV^'Hx,b = 0,p,fi) J (e^b.q _ (2.34) 

The Sudakov exponent due to perturbative evolution of the wave function, Eq. (p.5|), is determined 
to leading order in by the above expression. At this stage the running coupling is introduced. 
A perturbative analysis at the logarithmic level |^ indicates that the scale of the coupling is set by 
the transverse momentum of the emitted gluon and it encodes the information that the strength 
of the interactions increases for emission at large impact parameter. The Sudakov exponent now 
reads 



|q| Jh\ ) 



When the lower limit of the dq^ integration is set to 1/6^ the result in Eq. ( |2.6D is recovered. 
However, apart from the case of small size onia, for the usual light mesons 1/6 ~ Aqcd. It is then 



9 



apparent that after the introduction of the running coupHng the Sudakov exponent wih become 
dominated by the lower end of the dq^ integration where perturbation theory itself is ill defined. 
Nevertheless, what we are interested in here is not the numerical stability of the perturbative results 
but signals of nonperturbative corrections. We therefore proceed by allowing the infrared regulator 
to become ©(Aqcd) and introduce the Borel transformation of the Sudakov exponent defined as 

POO 

Swf{x,b,Q,fi;as) = da Swf{x,h,Q, n;a) e~''^"'' . (2.36) 
JO 

The Borel image Sy^f can be obtained by first using the following representation of the one loop 
running coupling: 

a,(q2)=/ dae-'^^oi'^^'i'/^'). (2.37) 
Jo 

Substituting the above expression into Eq. ( |2.35| ) and integrating over q we get 

Comparing the result with the definition ( p.36| ) we read off the Borel image 

~ Cf x{l-x)[p+f r(l - a/?o) 

^'"^^^^^ a/?or(l + a/3o)- ^'"^^^ 

The pole at o" = is of UV origin. Infrared sensitivity is parametrized by the presence of IR 
renormalon poles at cr = n/ (3q , n = 1, 2, .... The leading IR renormalon is at = l//3o which leads 
to 0(Aq(jj-)6^) power corrections of the form 

S^f = S^J + CA|cD&^ In ^^1-^ . (2.40) 

^^QCD 

Terms that are Q-independent have been omitted in the above equation. As in the case of the 
massive gluon scheme, it is this part of the coefficient that is n-independent to leading logarithmic 
in Q order and it arises from soft gluon emission, hence it is universal. The coefficient C contains the 
principal value integral over a, which is one possible way to define it. This ambiguity is anticipated 
since perturbation theory is ill defined and therefore it cannot predict the size of nonperturbative 



corrections, only their scaling with Q. We emphasize that the manipulations that led to Eq. ( 2.39 ) 
are formal and the transformation cannot be inverted back to as space without the use of some 
prescription (or convention). The value of the coefficient C will depend on this prescription. In 
addition, as with the massive gluon case, higher orders in ag are not suppressed, unless one invokes 
an additional assumption of say the freezing of the coupling at small momenta [T6| . 

The connection between the running coupling and the massive gluon calculation is estab- 
lished by comparing Eq. (pl0| ) with Eq. (|^). Note that in the massive gluon result we have not 
written the analogue of the S^J that we see above. This is because in the Mellin transformation 
analysis of the massive gluon integrals we looked only for contributions that are nonanalytic in 
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the regulator and not in the impact parameter 6^. Had we chosen the latter we would have 
generated the perturbative answer S^J in the massive gluon scheme as well. The nonperturbative 
pieces in the two calculations can be mapped onto each other, up to a multiplicative constant, via 
the identification 

A^lnA^^A^jCD- (2.41) 



In Eq. ( p. 31 ), the A^ In^ A^ term indicates that the coefficient of the term A^lnA^ will depend on 
InQ^, being the only other scale involved, hence we also identify (1/2)A^ In^ A^ ^ ^qcd ^^^qcd- 
The existence of A^ In A^ contributions is a signal for the presence of a condensate-like term. Had 
only this term been present then one could expect that the power corections could be captured from 
an O.P.E. in terms of local operators. The existence of A^ In^ A^ contributions, though, signals that 
such an expansion is not possible for the wave function of a light hadron. An operator expansion 
formalism for the wave function must necessarily involve nonlocal operators that extend along the 
light cone and, in analogy with the usual OPE, their expectations are to be parametrized by nonlocal 
condensates. The possibility of parametrizing power corrections to the meson electromagnetic form 
factor in terms of nonlocal condensates has been discussed in Ref. |17]. 



The other important piece of information in the leading power correction term is that its 
6^ dependence leads to a Gaussian factor for the meson wave function. This is obtained from 
Eqs. Q and ( pO] ). 

?PT „„„ / f ~ ^)Q ^\ a2 ;,2 

^QCD 



V = V^' exp [-Clni ^ ) ^^CD^' ) , (2-42) 



with 'P^'^ denoting the wave function perturbatively evolved with exponent given in Eq. ( |2.6D . We 
turn to a discussion of this form in the next section. 

3 The Gaussian form of the meson wave function 

The wave function in Eq. ( ^3T|) and Eq. (p^ ) can be rewritten in the following suggestive form: 



V{x, b, Q) = (/.(x, Q) exp (-5^J - h^S2{x{l - x)Q^) - b'^Si{xil - x)Q^) + ...) . (3.1) 

Infrared renormalons at integer values generate a expansion for the Sudakov exponent with Q 
dependent coefficients. This form is the same as the one seen in the Drell-Yan process at measured 



Q± and has been derived by Korchemsky and Sterman in Ref. [10|. These authors have also given 
an operator definition of the power corrections in terms of pairs of Wilson lines and their transverse 
derivatives. For the same Drell-Yan process, the Gaussian form in impact parameter space has been 
noted years ago and outside the renormalon context by Collins and Soper |l^] in their treatment 
of the infrared sensitivity of the Sudakov exponent. 
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The dimensionfull functions 5*2, 5*4, ... cannot be normalized within perturbation theory. 
However, we have obtained some information beyond the summation of logarithmic corrections. If 
we assume that in Eq. ( ^.42 ) the unknown coefficient C is a number of 0(1) and that our methods 
capture correctly the x-dependence of the leading nonpeturbative corrections then 

52(Q)~A^cdM^(1-^)Q'), (3.2) 

where only the universal piece is retained, as it has been argued in the previous section. For fixed 
X ^ 0, power corrections lead to increasing suppression of emission at large impact parameter with 
increasing Q. This suppression is in addition to the one generated by S^J . However, for fixed Q, 
it is seen that this additional suppression becomes weak in the end-point region. The end-point 
region is enhanced relative to the central region by the Gaussian dependence of the meson wave 
function. This is not surprising since it is in the end-point region that the effective hard scale 
x{l — x)Q'^ becomes small and power corrections become more important. We emphasize that all 
this is on top of the Sudakov suppression, 5^J, as mentioned above. In this respect, the pattern 
of power corrections is very similar to that observed, for example in event shape variables where 
the leading power corrections [^] also come from the end-point region (two jet limit in this case), 
which is itself Sudakov suppressed. 

For exclusive processes there exist in the literature various models for the wave function 
with which we may compare our results. In the case of the meson, the most popular is the one 
obtained from an oscillator model of two constituent quarks boosted in the light cone frame. This 
leads to a wave function of the form |19|, [20| 

P(x, 6, Q) = (j){x, Q) exp [s^J- <kl> x(l - x) 6^) , (3.3) 

where the average < k'j_ is the oscillator parameter. Note that this Gaussian form is such that 
it can interpolate between the perturbative tail of the wave function given by exp(— 5^J) and the 
nonperturbative region. The 0{b'^) term in the exponent has been found numerically important in 
order for Sudakov resummed perturbative expressions for the meson form factor to be applicable 



in the subasymptotic Q region |21|. It is therefore interesting to compare its 2;-dependence with 
the expression we get from perturbation theory. We observe that like in Eq. (|3.2|), the exponent 
in Eq. ( ^.3| ) enhances the end-point regions compared to the central region. However, from the 
renormalon based approach, which is rooted in perturbation theory, one can expect a dependence 
on x(l — x) which is logarithmic but not proportional to it. The analysis of nonperturbative 
corrections as signaled by perturbation theory and phenomenological models for the meson wave 
function in the nonperturbative regime have similar qualitative behavior although the end-point 
enhancement in the former case is much milder than in the latter. For our approach to give results 
in exact correspondence with, say, the oscillator model, it would require a resummation of the 
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higher order in Og contributions to the In term at the level of the exponent. These could turn 
the In x(l — x) into a. x^{l — x)" dependence. It is not obvious to us how such a resummation would 
be implemented. Recall that in the region where nonperturbative corrections are generated, the 
coupling is = 0(1) and there appears to be no small parameter around which to build and then 
resum a perturbative series. More importantly, it is not known how to define the massive gluon 
scheme for higher loop corrections, and for the running coupling scheme it is only an assumption 
that as = as{kj_) all the way into the deep infrared region. We can say for sure, though, that any 
answer obtained from perturbation theory would necessarily have Q dependence in the Gaussian 



form that it generates, unlike Eq. (3^) which is Q'hidependent beyond the logarithmic corrections 



residing in S^J. It should be pointed out here that, recently, renormalon based models have 
been used to predict the x-dependence of the higher twist structure functions in deeply inelastic 



scattering with some phenomenological success [ 23 1 and this was one of our motivations for pursuing 
the above analysis. 

4 Summary 

In this paper we studied the nonperturbative corrections to the meson wave function using the 
methods that have already been developed for semi-inclusive cases, namely the massive gluon and 
the renormalon methods. We found that the leading nonperturbative corrections are of order 
0(Aqqj-,6^) for the two-quark wave function at transverse separation b. The exponentiation of 
such contributions leads to a Gaussian factor in addition to the Sudakov resummed logarithmic 
enhancements. Of particular interest is the x-dependence of this Gaussian factor. It leads to the 
conclusion that the power corrections arise from the end-point regions. This x-dependence has 
been compared with low energy Gaussian models for the wave function. 

It must be emphasized that both methods for obtaining the leading nonperturbative cor- 
rections have their origin in perturbation theory. They predict correctly the type of the power 
correction but since they are applied in a region where the coupling is normalized at low scales they 
have limited predictability for the coefficient of the power corrections unless additional assump- 
tions are introduced, such as freezing of the coupling. Thus our predictions for the x-dependence 
should be considered as another model. Similar approach for the x-dependence of the higher twist 
structure functions in D.I.S. has met with some phenomenological success. It would be interesting 
to apply this model to the phenomenological study of the meson electromagnetic form factor and 
elastic scattering. 
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